Experimental Constraints on the Spin and Parity of the A c (2880) 
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We report the results of several studies of the A+7r + 7r~X final state in continuum e + e~ annihi- 
lation data collected by the Belle detector. An analysis of angular distributions in A c (2880) + — > 
E c (2455)°' ++ 7T + ' _ decays strongly favors a A c (2880) + spin assignment of | over | or |. We find 
evidence for A c (2880) + -» E c (2520) 0,++ 7r + '~ decay and measure the ratio of A c (2880) + partial 
widths r(s°(2455)^) = 0.225 ±0.062 ±0.025. This value favors the A c (2880) + spin-parity assignment 
of | + over |~. We also report the first observation of A c (2940) + — > E c (2455)°' ++ 7r + '~ decay and 
measure A c (2880) + and A c (2940) + parameters. These studies are based on a 553 fb -1 data sample 
collected at or near the T(4S') resonance, at the KEKB collider. 

PACS numbers: 13.30.Eg, 14.20.Lq 



Charmed baryon spectroscopy provides an excellent 
laboratory to study the dynamics of a light diquark in 
the environment of a heavy quark, allowing the predic- 
tions of different theoretical approaches to be tested 0, 
0) B El- There are twelve experimentally observed 
charmed baryons for which the spins and parities are as- 
signed [f| They include ground states, spin excita- 
tions and lowest orbital excitations. Except for the A+, 
the J p quantum numbers for these states have not been 
determined experimentally but are assigned based on the 
quark model predictions for their masses. There are also 
six charmed baryons, recently observed at the CLEO 0], 
Belle [1, [t| and BaBar [l(| experiments, for which the 
spins and parities are not well constrained. The new 
states are in a mass region where the quark model pre- 
dicts many levels with small spacing, which makes the 
J p assignment difficult. In this Letter we investigate 
possible spin and p arity values of one such state, the 
A c (2880) + baryon [7|,ll0(, by studying the resonant struc- 
ture of A c (2880) + — > A+7T+7T - decays and performing an 
angular analysis of A c (2880)+ -> £ c (2455)°' ++ 7r+'- de- 
cays. We also report the first observation of A c (2940) + — > 
S c (2455)°' ++ 7r+'- decay and measure A c (2880)+ and 
A c (2940) + parameters. 

We use a 553 fb _1 data sample collected with the Belle 
detector at or 60MeV below the T(4S) resonance, at the 
KEKB asymmetric-energy (3.5 GeV on 8.0 GeV) e + e~ 
collider [ll[ . The Belle detector [l^] is a large-solid-angle 
magnetic spectrometer that consists of a silicon vertex 
detector, a 50-layer cylindrical drift chamber, an array 
of aerogel threshold Cherenkov counters, a barrel- like 
array of time-of-fiight scintillation counters, and an ar- 



ray of CsI(Tl) crystals located inside a superconducting 
solenoidal coil that produces a 1.5T magnetic field. An 
iron flux return located outside the coil is instrumented 
to detect muons and mesons. We use a GEANT 
based Monte-Carlo (MC) simulation [l3[ to model the 
response of the detector and to determine its acceptance. 
Signal MC events are generated with experimental run 
dependence in proportion to the relative luminosities of 
different running periods. 

A+ baryons are reconstructed using the pK~n + decay 
mode (the inclusion of charge conjugate modes is im- 
plied throughout this Letter). To select proton, charged 
kaon and pion candidates we use the same track qual- 
ity and particle identification criteria as for observa- 
tion of the E c (2800) isotriplet [||. The invariant mass 
of the pK~n + combinations is required to be within 
±8MeV/c 2 (1.6a) of the A+ mass value, recently mea- 
sured by BaBar 14] . To improve the accuracy of the A+ 
momentum measurement we perform a mass constrained 
fit to the pK~ir + vertex. We combine A+ candidates 
with the remaining tt + tt~ candidates in the event. To re- 
duce the combinatorial background we impose a require- 
ment on the scaled momentum of the A+ 7r + 7r~ combina- 
tion x p = P*/\/E^ am - M 2 > 0.7, where p* is the mo- 
mentum and M is the invariant mass of the combination, 
-^beam ^ s the beam energy, all variables being measured in 
the center-of-mass frame. To improve the M(A+7r + 7r~) 
resolution we perform an interaction point constrained 
fit to the A+7r + 7r~ vertex. 

To measure the A c (2880) + mass and width we ap- 
ply an additional requirement that cither M (A+7T - ) or 
M(A+7r + ) be in the S c (2455) signal region defined as 
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2450MeV/c 2 < M < 2458MeV/c 2 . Whereas 35% of sig- 
nal events pass this cut, only 12% of background events 
do so. From MC simulation we find that the mass res- 
olution for the A c (2880)+ -> E c (2455)°' + +7r+<- decays 
depends strongly on the decay angle 9, defined as the an- 
gle between the pion momentum in the A c (2880) + rest 
frame and the boost direction of the A c (2880) + . To as- 
sure good resolution for the A c (2880) + mass and width 
measurement we require cos 9 > 0. This requirement also 
helps to suppress combinatorial background. The result- 
ing M(A+7r + 7r - ) distribution is shown in Fig. [1] One 




M(A+7r+7T-), GeV/c 2 

FIG. 1: The invariant mass of the A^tt + tv~ combinations for 
the E c (2455) signal region (histogram) and scaled sidebands 
(dots with error bars). The fit result (solid curve) and its 
combinatorial component (dashed curve) are also presented. 

can see clear peaks from the A c (2765) + and A c (2880) + . 
A peak in the region M = 2940MeV/c 2 is associated 
with the A c (2940) + baryon recently observed in the D°p 
final state by BaBar [lOj. Scaled S c (2455) sidebands, 
which are also shown in Fig. [TJ are featureless in the 
region of the A c (2940)+. The £ c (2455) sidebands are 
defined as 2438MeV/c 2 < Af(A+7r) < 2446MeV/c 2 and 
2462MeV/c 2 < M(A+tt) < 2470MeV/c 2 . 

We perform a binned likelihood fit to the A+7r + 7r~ 
mass spectrum of Fig. Q] to extract the parameters and 
yields of the A c (2880)+ and A c (2940)+. The fitting func- 
tion is a sum of three components: A c (2880) + signal, 
A c (2940) + signal and combinatorial background func- 
tions. As shown below, the favored spin-parity assign- 
ment for the A c (2880)+ is | + , therefore the A c (2880)+ 
signal is parameterized by an F-wave Breit-Wigner func- 
tion convolved with the detector resolution function, de- 
termined from MC (a = 2.2MeV/c 2 ). The A c (2940)+ 
signal is an S-wave Breit-Wigner function convolved with 
the detector resolution function (a = 2.4MeV/c 2 ). The 
background is parameterized by a third-order polyno- 
mial. The fit is shown in Fig. [IJ and the results are sum- 
marized in Table |TJ The signal yield is defined as the in- 
tegral of the Breit-Wigner function over a ±2.5r interval. 



TABLE I: Signal yield, mass and width for the A c (2880) + and 
Ac (2940) + . The first uncertainty is statistical, the second one 
systematic. 



State 


Yield 


M, MeV/c 2 


T, MeV 


A c (2880)+ 
A c (2940)+ 


690 ± 50 
220± 8 6 ° 


2881.2 ± 0.2 ± 0.4 
2938.0 ± 1.3±l ° 


5.8 ±0.7 ±1.1 
13tgt a ? 



The normalized \ 2 of the fit is x 2 /d.o./. = 132.2/134. If 
the A c (2940) + signal is removed from the fit, the double 
log likelihood changes by 59.8, which corresponds (for 3 
degrees of freedom) to a signal significance of 7.2 stan- 
dard deviations. 

To estimate the systematic uncertainty on the re- 
sults of the fit we vary the background parameteriza- 
tion, using a fourth-order polynomial and the inverse 
of a third-order polynomial. We include the A c (2765) + 
signal region into the fit interval, parameterizing the 
A c (2765) + signal by an S-wave Breit-Wigner function. 
The A c (2765) + mass and width determined from the fit 
are M = (2761 ± 1) MeV/c 2 and T = (73 ± 5) MeV. We 
vary the selection requirements; we take into account the 
uncertainty in the A+ mass of ±0.14MeV/c 2 [lj], the 
mass scale uncertainty of 1q]21 MeV/c 2 [IH and the un- 
certainty in the detector resolution of ±10% as estimated 
by comparison of the inclusive A+ — > pK~ir + signal in 
data and MC. In the region between the A c (2880) + and 
A c (2940) + signals the fit is systematically below the data 
points, which might be due to a presence of an additional 
resonance or due to interference. We take into account 
these possibilities as a systematic uncertainty. In each 
case we consider the largest positive and negative varia- 
tion in the A c (2880)+ and A c (2940) + parameters to be 
the systematic uncertainty from this source; each term 
is then added in quadrature to give the total system- 
atic uncertainty, quoted in Table HI The main sources 
of the systematic uncertainty are a possible contribution 
of the A c (2765) + tail into the fit region (the shape of 
the tail is not well constrained) and the excess of events 
between the A c (2880)+ and A c (2940)+ signals. None of 
the variations in the analysis alters the A c (2940) + signal 
significance to less than 6.2 standard deviations. 

For further analysis, we remove the cos 9 > re- 
quirement. To study the resonant structure of the 
A c (2880) + — > A+7T+7T - decays we fit the A+7r+7r~ mass 
spectrum in Af (A+7T ) bins. By isospin symmetry, we ex- 
pect equally many decays to proceed via a doubly charged 
E c (2455) (E c (2520)) as via a neutral one. Since the 
corresponding doubly charged and neutral channels are 
kinematically separated in phase space, we combine the 
Af (A+7T + 7T~ ) distributions for Af (A+tt - ) and Af(A+7r+) 
bins. To fit the A+7r + 7r~ mass spectra we use the same 
fit function as described above. The A c (2880) + and 
A c (2940) + parameters are fixed to the values in Tablefl] 
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The A c (2880)+ yield as a function of M(A+tt ± ) is shown 
in Fig. [2j We find a clear signal for the X c (2455) and 
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FIG. 2: The A c (2880) + yield as a function of M(A+tt ± ). The 
histogram represents the result of the fit. 

an excess of events in the region of the E c (2520). We 
perform a \ 2 fit to the A+7r ± mass spectrum of Fig. [2] to 
extract the yields of the £ c (2455) and S c (2520). The fit- 
ting function is a sum of three components: £ c (2455) 
signal, £ c (2520) signal and a non-resonant contribu- 
tion. The £ c (2455) and £ c (2520) signals are parameter- 
ized by a P-wave Breit-Wigner function convolved with 
the detector resolution functions, determined from MC 
(a = 0.9MeV/c 2 for the £ c (2455) and a = 1.5MeV/c 2 
for the £ c (2520)). The mass and width of the £ c (2455) 
are floated, while the mass and width of the £ c (2520) are 
fixed to the world average values The shape of the 
non-resonant contribution is determined from MC assum- 
ing a uniform distribution of the signal over phase space. 
The fit is shown in Fig. [2] We find the ratios of A c (2880)+ 
partial widths E^^^f] _ Q 4Q4 ± Q Q21 ± Qm ^ 

r(s (2520)^) = .091±0.025±0.010and F ^f™K\ = 

r(A+7r+7r-) r(S c (2455)7r±) 

0.225 ± 0.062 ± 0.025, where the uncertainties are statis- 
tical and systematic, respectively. The S c (2455) param- 
eters determined from the fit M = (2453.7±0.1) MeV/c 2 
and r = (2.0 ± 0.2) MeV are consistent with the world 
average values The normalized x 2 °f the fit is 

X 2 /d.o.f. = 106.6/75. The significance of the S c (2520) 
signal is 3.7 standard deviations. 

To estimate the systematic uncertainties on the ratios 
of A c (2880) + partial widths we vary the A c (2880)+ pa- 
rameters, fit interval and background parameterization 
in the fit to the M(A+7r + 7r~) spectrum; we vary the 
E c (2520) parameters; we allow the shape of the non- 
resonant contribution to float in the fit, parameterizing 
it with a second-order polynomial multiplied by a thresh- 
old function or by a third-order polynomial; we take into 
account the uncertainty in the detector resolution and in 
the reconstruction efficiency. None of the variations re- 



duces the significance of the S c (2520) signal below three 
standard deviations. 

To perform angular analysis of A c (2880) + — > 
£ c (2455) 0,++ 7r + '~ decays we fit the A+7r + 7r~ spectrum 
in cos# and (f> bins for the S c (2455) signal region and 
sidebands. Here, 4> is the angle between the e + e~ — > 
A c (2880) + A reaction plane and the plane defined by the 
pion momentum and the A c (2880) + boost direction in 
the rest frame of the A c (2880) + . Figure [3] shows the 
yield of A c (2880) + as a function of cos8 and <j>, after 
S c (2455) sideband subtraction (to account for nonreso- 
nant A+7r + 7r~ decays) and efficiency correction. 
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FIG. 3: The yield of A c (2880)+ -» E c (2455)°' ++ 7r+ " decays 
as a function of cos 8 and <t>. The fits are described in the text. 



The parameterization of A c (2880) + — > E c (2455)7r de- 
cay angular distributions depends on the spin of the 
A c (2880) + . For the spin | hypothesis both cos# and 
<j) distributions are expected to be uniform [l6l |. \ 2 fits 
to a constant are shown in Fig. [3] by a dotted line. The 
agreement is good for tfi: x 2 /d.o.f. = 5.3/9, but poor for 
cos 9: x 2 /d-o.f. = 46.7/9. 

The angular distribution for the spin | hypothesis 
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is hf 



IP 



3/2 



3 1 
— [p 3 3 sin 2 + pn(- +cos 2 i 

47T O 



2,2 

— =Rep 3 _i sin 9 cos 20 -=Rep 3 i sin 29 cos <p\ , 

V 3 V 3 

where pij are the elements of the production density 
matrix. The diagonal elements are real and satisfy 
2(/5 3 3 + pix) = 1. Since the measured distribution in 
<j) is consistent with being uniform (this also holds sep- 
arately for cos 6 > and cos^ < samples), the non- 
diagonal elements are small. The result of the fit to 
the cos 9 spectrum for the spin | hypothesis is shown 
in Fig. [3] with a dashed curve. The agreement is poor: 
X 2 /d.o.f. = 35.1/8. 

The angular distribution for the spin | hypothesis 
is M 



W 5/2 = -[p552(5cos 4 6»-2cos 2 6» + l) + 
8 

p 33 (-15cos 4 9 + 14cos 2 9 + 1) + p u 5(l - cos 2 9)% 

where non-diagonal elements are ignored. The result of 
the fit to the cos6> spectrum for the spin | hypothesis 
is shown in Fig. [3] with a solid curve. The agreement is 
good: x 2 /d-o.f. = 12.1/7. We find p 55 = 0.09±0.02 and 
p 33 = 0.00 ±0.03. Thus the A c (2880)+ populates mainly 
thehelicityzb 1 ; states, 2p xl = l-2p 33 -2p 55 = 0.82±0.05. 

The x 2 difference of the spin i (|) and spin | fits is dis- 
tributed as x 2 with two degrees (one degree) of freedom, 
therefore the exclusion level of the spin i (|) hypothesis 
is 5.5 (4.8) standard deviations. 

To estimate the systematic uncertainty in the angular 



analysis of the A c (2880)" 



S c (2455) 



0,++^+,- 



decay 



we vary the A c (2880) + parameters, fit interval and back- 
ground parameterization in the fit to the M(A 



+ 7T + 7T 



_ ) 

spectrum. None of the variations alters the exclusion 
level of the spin i (|) hypothesis to less than 5.5 (4.5) 
standard deviations. 

The Capstick-Isgur quark model predicts the lowest 
J p = |~ A+ state at 2900MeV/c 2 and the lowest J p = 
| + A+ state at 2910MeV/c 2 Q. The typical accuracy 
of quark model predictions is 50MeV/c 2 , therefore the 
agreement with the experimental value for the A c (2880) + 
mass is quite good. The lowest spin | states are well 
separated from the next J — | levels (3130MeV/c 2 for 
negative and 3140MeV/c 2 for positive parities) and from 
J = f levels (3125 MeV/c 2 for negative and 3175 MeV/c 2 
for positive parities). 

Heavy Quark Symmetry predicts R = 



1.4 

5 + 



r(£ c (2520)7r) _ 
r(S c (2455)ir) ~~ 

0.36 for the 



for the | state and R = 0.23 
state @, E3]. The measured value R = 0.225 ± 



0.062 ± 0.025 favors the positive parity assignment for 
the A c (2880)+. 



The | + assignment for the A c (2880) + makes it a 
special state that lies on the leading A+ Regge trajec- 
tory, whose lower J p members are the | A, 



and | 



A c (2625) + . The § assignment for the A c (2880)+ based 
on a string model for baryons was proposed in Ref. [l8j |. 

In summary, from angular analysis of A c (2880) + — > 
S c (2455)°' ++ 7r+-- decays we find that a A c (2880)+ spin 
hypothesis of | is strongly favored over i and |. We 
find first evidence for S c (2520)7r intermediate states 
in the A c (2880) + — > A+7T+7T - decays and measure 

r(S°(2455)^±) = °' 225 ± °- 062 ± °- 025 ' This value is 

in agreement with Heavy Quark Symmetry predictions 



and favors the | 



over the | hypothesis for the spin- 



parity of the A c (2880) + . We also report the first obser- 
vation of A c (2940) + — * E c (2455)7r decays, and measure 
the A c (2880)+ and A c (2940)+ parameters. 
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